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Understanding  the  Effect  of  Atmospheric 
Turbulence  on  Optical  and  Infrared  Propagation 
using  Hilbert  Phase  Analysis 

C.  Font/  G.C.  Gilbreath/  M.P.J.L.  Chang/  and 
E.S.  Oh^ 

^Information  Technology  Division 
^University  of  Puerto  Rico,  Mayaguez  Campus 
^Remote  Sensing  Division 

Introduction:  The  refractive  index  structure 
parameter  is  a  figure  of  merit  used  to  describe  the 
magnitude  of  the  effect  of  turbulence  in  the  atmo¬ 
sphere.  Understanding  the  influence  of  different  local 
climate  parameters  upon  is  an  important  step 
toward  developing  a  prediction  model  for  in-field  laser 
interrogators.  The  layered  structure  of  the  atmosphere 
and  its  random  behavior  present  unique  challenges  to 
theoretical  and  experimental  studies  in  propagation 
research.  This  article  presents  a  novel  technique  termed 
Hilbert  Phase  Analysis  (HPA)  to  provide  insights  into 
the  physical  interactions.  We  explore  the  technique’s 
use  in  analyzing  the  dependence  of  on  local  climate 
variables  using  data  recently  obtained  from  campaigns 
conducted  in  Puerto  Rico. 

The  Method:  HPA  is  a  new  technique  that  consists 
of  Empirical  Mode  Decomposition  (EMD)  and  the 
Hilbert  Transform  (HT).  EMD  is  an  instantaneous 
frequency  hlter  bank  that  generates  a  family  of  empiri¬ 
cal  eigenmodes  called  Intrinsic  Mode  Eunctions  (IMF) 
from  a  non-linear,  non-stationary  time  series.^  By 
transforming  the  IMFs  with  HT,  we  produce  physically 
interpretable  instantaneous  phase  angles,  tp,-.  These 
values  are  akin  to  Fourier  phase  angles,  except  that 
they  apply  per  sample  point.  This  feature  now  enables 


the  study  of  the  temporal  phase  relationship  for  many 
measurables  in  a  point-wise  manner.  This  is  not  pos¬ 
sible  using  standard  Fourier  techniques.  Legacy  time 
series  tools  such  as  wavelets  provide  relatively  poor 
temporal  resolution.  Figure  8  illustrates  how,  by  using 
HPA,  the  phases  (pi  and  (p2  of  two  interdependent 
signals  can  be  shown  to  be  related.  On  the  left.  Signal 
1  ’s  vector  rotates  clockwise,  while  that  of  Signal  2  on 
the  right  rotates  counter-clockwise.  The  time-depen¬ 
dent  phase  angles  of  these  signals  are  clearly  inversely 
correlated.  The  relationship  can  provide  insight  into 
how  climatic  parameters  might  affect  C^. 

Experiment:  A  series  of  campaigns  to  collect 
and  simultaneous  climate  data  are  being  performed  at 
the  University  of  Puerto  Rico  at  Mayagiiez.^  The  results 
described  here  are  from  a  campaign  conducted  in 
February  and  March  2006.  Scintillometers  were  placed 
on  the  roofs  of  buildings  at  the  main  campus  of  the 
university.  The  transmitter  and  receiver  were  separated 
by  90  meters  over  land.  Climate  data  was  acquired  at  a 
nearby  weather  station. 

Data  Analysis:  Figure  9(a)  compares  meteoro¬ 
logical  data  with  over  time.  This  data  is  from  a 
very  clear  and  sunny  day.  We  can  see  close  tracking 
between  the  measured  solar  radiation  flux  (SR)  and 
the  temperature  (T),  which  is  expected,  since  the  solar 
flux  controls  the  heating  of  the  atmosphere.  We  also 
observe  similar  tracking  of  with  T,  validating  that 
temperature  plays  a  major  role  in  defining  optical  tur¬ 
bulence.  From  the  relative  humidity  (RH)  record,  we 
see  an  inverted  relationship  to  T  and  therefore  to  . 

If  we  directly  cross-correlate  and  a  climate  variable, 
the  resulting  function  is  very  difficult  to  interpret,  since 
Cn  is  highly  variable  with  time. 
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FIGURE  8 

Illustration  of  HPT  phase  angles 
where  the  colors  represent  the 
changes  in  time.  The  arrows 
show  the  signal  vectors,  which 
change  over  20  time  steps. 
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Solar  Radiation  vs.  Time  x10“^^ 


Cross  Correlation  of  HPA 


Lag  (10  mins/interval) 

(b) 

FIGURE  9 

(a)  Data  for  March  9,  2006;  (b)  Cross-correlations  of  the  instantaneous  phase  angles  of  the  different  measurements 
within  ±300  minutes  of  the  zero  lag  position.  Black:  The  correlation  function  between  solar  radiation  flux  (SR)  and 
temperature  (T)  phase  angles,  showing  a  strong  positive  correlation.  Red:  The  correlation  function  between  Cp  and  tem¬ 
perature  (T)  also  showing  a  strong  positive  correlation.  Blue:  The  correlation  function  between  €„  and  relative  humidity 
(RH),  showing  a  negative  correlation. 
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Applying  HPA  to  this  family  of  measurements 
provides  us  with  the  time-varying  phase  angles.  Since 
we  expect  correlations  between  our  raw  parameter 
records,  we  calculated  cross-correlations  of  the  HPA- 
determined  phase  angles,  as  shown  in  Fig.  9(b).  The 
curve  plotted  in  black  is  the  cross-correlation  function 
of  the  instantaneous  phase  angles  SR  and  T;  it  shows 
that  T  faithfully  follows  SR,  with  an  approximate 
ten-minute  lag.  The  cross-correlation  function  of 
and  T  indicates  that  the  turbulence  parameter  follows  a 
similar  relationship.  These  results  support  the  conven¬ 
tional  theory  that  temperature  strongly  influences  , 
thereby  calibrating  the  HPA  technique  against  known 
atmospheric  behaviors.  Finally,  the  cross-correlation 
function  between  and  RH  shows  an  inverse  depen¬ 
dence,  which  corroborates  the  first-order  approxima¬ 
tion  of  an  inverse  relationship  existing  between  tem¬ 
perature  and  relative  humidity. 

Conclusions:  HPA  provides  the  ability  to  study 
the  time-varying  phase  angles  of  non-linear  signals. 


In  this  research,  the  technique  was  tested  with  optical 
turbulence  data  and  its  simultaneous  local  climate 
record  taken  in  a  recent  campaign  in  Puerto  Rico.  The 
technique  demonstrated  consistency  with  conventional 
physics.  That  is,  the  direct  influence  of  temperature  on 
Cn  was  verified.  Additionally,  the  technique  reveals 
an  approximate  inverse  relationship  between  and 
relative  humidity.  The  experiment  is  being  extended 
to  over-the-water  ranges  of  up  to  600  meters  between 
Magueyes  Island  and  the  Villa  Parguera  Hotel  in  La 
Parguera,  Lajas,  Puerto  Rico,  as  shown  in  Fig.  10. 
[Sponsored  by  ONR] 
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FIGURE  10 

Propagation  path  shown  for 
Cn  experimental  studies  for 
2006-2007  campaigns. 
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